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The photon spectrum in orthopositronium — > 37 decay is computed using effective field theory 
methods. For energetic photons, the spectrum agrees with the Ore-Powell result, but deviates from 
it when the photon energy is comparable to the positronium binding energy. The decay spectrum in 
this region depends on a positronium structure function, which is computed in this paper. At still 
lower energies the photon spectrum is dominated by the parapositronium resonance contribution. 
Our results are compatible with Low's theorem on soft photon emission. 
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I. INTRODUCTION 

Since its discovery in 1951 pj Positronium has be- 
come a high precision laboratory for testing QED. A de- 
tailed study of the properties of the e + e~ spin triplet 
(orthopositronium: o-Ps) and spin singlet (parapositro- 
nium: p-Ps) states has provided impressive confirmation 
of QED radiative corrections. 

Recently, it was pointed out |2j that the calculation 
of the o-Ps decay spectrum, first performed by Ore and 
Powell 3], is in apparent contradiction with Low's theo- 
rem, i.e. that the soft photon spectrum is not consistent 
with QED gauge invariance in combination with analyt- 
icity. While Low's theorem applied to this decay requires 
that the spectrum vanishes as £|2, with E 1 the energy of 
the radiated photon, the standard Ore-Powell calcula- 
tion predicts a 0{E 1 ) behavior for E 1 — > 0. According 
to the authors of Ref. [2|, this failure is traced back to 
the different radiation properties of the free charged lep- 
tons which are used as the asymptotic states in typical 
o-Ps— ► 37 calculations, and those of positronium which 
is a neutral spin-1 boson. 

In this paper we compute the orthopositronium — ► 37 
decay spectrum using non-relativistic effective field the- 
ory (NREFT) methods developed to study bound states 
in QED and QCD 0,HIE0- The ex P ansion parameter 
of the effective theory is the velocity v <C 1 of the electron 
in the Ps bound state. [We will use the velocity power 
counting of Refs. and treat a ~ v.] We show that 

the Ore-Powell calculation can be reconciled with Low's 
theorem if the latter is carefully applied considering all 
the energy scales present in the problem. 

The three important scales for the photon decay spec- 
trum are the electron mass to, the binding energy of order 
ma , and the hyperfine splitting between the singlet and 
triplet states of order ma 4 . The Ore-Powell computation 
is valid for photon energies E 1 ^> ma 2 . When E 1 is of 



order 



the Ps binding energy can not be neglected. 



The decay amplitude depends on a sum over an infinite 
set of excited Ps states, and can be written in terms of a 
Ps structure function. We will compute the Ps structure 
function using NRQED to leading order in the v expan- 
sion. When E 1 is of order ma 4 , the decay amplitude is 



dominated by the o-Ps — > p-Ps transition. We will show 
that the photon spectrum crosses over from an E^ behav- 
ior above the structure function region to an E^ behavior 
below the p-Ps resonance region. 

The outline of the paper is as follows. In section [D] 
we review Low's theorem for soft photon emission and 
its implications for Ps decay. In section ITTTl we summa- 
rize the standard Ore-Powell computation of the decay 
spectrum. The NREFT calculation is discussed in the 
next three sections. Section IIVI discusses the matching 
between QED and NRQED, section IV! computes the de- 
cay amplitude in NRQED to leading order in v in terms 
of Ps matrix elements, and section IVTl discusses the com- 
putation of the required matrix elements. Results and 
conclusions are given in Section lVllI 



II. LOW'S THEOREM AND O-PS -» 3 7 DECAY 

Low's theorem [j| gives the amplitude for soft photon 
emission in the scattering of charged particles. It states 
that the first two terms of the series expansion in powers 
of the photon energy of a radiative amplitude X — > Y7 
may be obtained from a knowledge of the corresponding 
nonradiative amplitude X — > Y. Expanding the radia- 
tive amplitude AP in the limit k — > 0, 



Mo 



+ M 1 + 0(k), 



(1) 



where e is the photon polarization, one finds that Afo 
and Mi are independent of k and completely determined 
from the nonradiative amplitude To, its derivatives in 
physically allowed kinematic directions, and the electro- 
magnetic properties of the particles involved Q- The 
terms A-^0,1 are given by pole diagrams. The 0(k) am- 
plitude has both pole and non-pole contributions. 

The term A^o arises from the emission of a photon by 
ingoing or outgoing charged particles and is proportional 
to To times the universal factor —Qi e- Pi/k ■ pi, summed 
for all the external lines in the diagram. Note that if the 
nonradiative process involves no moving charged parti- 
cles or is forbidden due to some selection rule, then Mo 
is identically zero. The term Mi can be expressed as a 



2 



function of the magnetic moments, the amplitude To and 
its derivatives with respect to internal variables which are 
not subject to constraint, e.g. the energy and angles. The 
crucial observation is that unphysical derivatives with re- 
spect to masses (pf) can be shown to cancel out 0. 

Combining the amplitude behavior with that of the 
phase-space, the low-frequency form of the photon spec- 
trum is 



Pi 



dr 



A_ 

E~, 



B + G(T 7 ), 



(2) 



where A is proportional to \M.q\ and B is the MqM\ 
interference term. If Mo vanishes, the soft photon decay 
spectrum is of order T 7 di? 7 ; if both A4q and A'Jivanish, 
it is of order E^ dE y . 

In the three-photon decay of o-Ps, one of the photons 
can have an arbitrarily small energy. The process can 
be then viewed as the radiative version of the o-Ps^ 2j 
decay. As the two-photon decay of o-Ps is not allowed 
by charge conjugation invariance, the direct application 
of Low's theorem yields A^o.i = so that the o-Ps^ 37 
amplitude is of order 0(E 7 ), and the decay spectrum is 



dr o 



Ps->37 



d£L 



Ei 



(3) 



as E-y 



0. This is in contradiction with the Ore-Powell 



spectrum (see Eq. (11 81) ^) which vanishes linearly with E 1 , 
and is the contradiction pointed out in Ref. [2| • 

To understand the origin of the contradiction, it is 
worth noting that in the derivation of Low's theorem, 
one takes the limit E 7 — > and neglects all states other 
than those degenerate with the incoming and outgoing 
states, i.e. one uses E 1 <C AE, where AE is the en- 
ergy gap to excited states. The amplitudes Mo,i depend 
on the charge and magnetic moment couplings between 
all states degenerate with the initial or final states. A 
more general version of Low's theorem gives the decay 
spectrum for small E-y without taking the strict E 7 — > 
limit. One treats all states with AE -C E 1 as degenerate 
states, and includes them in the computation of charge 
and magnetic moment matrix elements for the purposes 
of Low's result. Which states are included in Low's the- 
orem then depends on the magnitude of E 7 . 

In the case of Ps decay, consider the case where the 
photon energy is much larger than the binding energy. 
Then all Ps states (including o-Ps, p-Ps, radial excita- 
tions, etc.) are degenerate for the purposes of Low's the- 
orem. In this case, there is a non-zero magnetic dipole 
matrix element between o-Ps and p-Ps, so that Aii does 
not vanish in this extended space of states. As a result, 
the decay spectrum vanishes linearly with E 7 . This is the 
approximation under which the Ore-Powell calculation is 
valid, as we will see in more detail later. For energies 
much smaller than the o-Ps-p-Ps hyperfine splitting, p- 
Ps as well as radial excitations are treated as excited 
states, the matrix element A^i vanishes, and the spec- 
trum is of order E^,. 
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FIG. 1: Three photon annihilation graph. The graph is 
summed over the 3! permutations of the photons. 



III. ORTHOPOSITRONIUM DECAY 
AMPLITUDE 

The Ore-Powell calculation Q of the 37 annihilation 
of o-Ps is given in many textbooks. We briefly review 
the derivation here, using the notation of Ref. LLfJ- The 
amplitude for an electron and positron to annihilate into 
three photons is given, to lowest order in a, by the graph 
in Fig.^ In the center-of-mass frame, the momentum of 
the electron is p and that of the positron is —p. The on- 
shell decay amplitude is a function of p. In the NRQED 
power counting p is of order 17 so one can expand the 
amplitude in a power series in v. The leading term in 
the expansion, which is all that is required here, is the 
p — * limit of the annihilation amplitude. 

Pick a gauge where the photon polarizations are 
h = 0, e? = 0. The sum of all 



purely transverse, Q • 
six QED graphs gives 



.4 
V 



V 



2m 2 ' 

(e 2 • e 3 ) ei + (e x • e 3 ) e 2 + (ei • e 2 ) e 3 
+ (6 2 ■ e 3 ) di + (St ■ e 3 ) 8 2 - (<5i • S 2 ) e 3 
+ (S 3 ■ ex) 8 2 + (8 2 ■ e%) 8 3 - (8 2 ■ <5 3 ) ei 
+ (<5i • e 2 ) 8 3 + (8 3 ■ e 2 ) £1 - (<5 3 • <5i) e 2 , 



where 



Si 



(4) 



(5) 



and we use non relativistic Pauli spinors 4> , x for the 
electron and positron respectively, 







(6) 



Let us take fc 3 as the photon with vanishing energy. As 
E 3 — > 0, the high-energy photons become antiparallel, so 
we have fci = — fc 2 . In this limit, Eq. |0J becomes 

V = 2 (e 2 • e 3 ) ei + 2 (ei • e 3 ) e 2 



3 



+2 



2 (ei • e 2 ) (ki ■ e 3 ) ki 
2 [ffei • £3) fe 2 • fc 3 J - f fei • fe 3 J (e 2 • e 3 
2 ( fc 3 • fci) (e 3 • e x ) - (k 3 ■ ei^J (e 3 ■ hi 
ei • fe 3 ) (e 2 • e 3 ) - (ei ■ e 3 ) (e 2 • fc 3 



£2 
fei, 
(7) 



since ei, 2 • fci = 0, and we have used the identity 
[(a x b) ■ c] [(sc xy)>z] 
= (a • sc) (6 • y) (c • z) + (a • y) (b ■ z)(c- x) 
+ (a ■ z)(b- x) (c • y) - (a ■ z) (b ■ y) (c ■ x) 
- (a • y) (6 • a;) (c ■ z) - (a ■ x) (6 ■ z) (c ■ y) 



A. Ore-Powell Decay Spectrum 



(8) 



2.5 




The o-Ps annihilation amplitude can be obtained from 
the free-particle decay amplitude in Eq. (0J by taking 
the matrix element between the o-Ps state and the vac- 
uum. If the o-Ps momentum space wavefunction is 4> (p) 
and the annihilation amplitude is ^4(p), the bound-state 
decay amplitude is 

^A( P )MP)- 0) 

To lowest order in v, A(p) is a constant, A(0), and Eq. © 
becomes 

A(0)M0), (io) 

where ip o {0) is the o-Ps position space wavefunction at 
the origin, x = 0. Using 

1 _,„,/„ 2 



-0o(x) 



-|x|/o 



(11) 



for the IS wavefunction, spin-averaging \A\ 2 in Eq. Q 
(see Ref. 0|) gives, 

JEW 4e ° 



spin 



3m z 



|^ o (0)| 2 ^ (l-cos0 12 ) 2 , (12) 

cyclic 



where cos#i 2 = fci • fe 2 . In terms of the dimensionless 
variables = Ei/m, 

x 1 + x 3 — 1 



(1 — cos ( 



xia; 3 



(13) 



and the differential decay rate is 

[xi +x 3 - l) 2 



dr 37 = n _ dzidcc 3 



9tt 



'1*3 



+ 



(l-^l) 2 



+ 



(1-Z 3 ) 2 



x 3 (2 — — x 3 ) 2 x{(2 — xi—x 3 ) 



(14) 



FIG. 2: The Ore-Powell orthopositronium decay spectrum. 



using the phase space factor 
1 2(2m) 2 



6 256tt 3 



&x\dx 3 . 



(15) 



In addition to the usual three-body phase space, we have 
1/6 for three identical particles in the final state, and 
2(2m) to convert from the non-relativistic normalization 
of the o-Ps state to the relativistic normalization. 

For the decay spectrum we can perform the integration 
over x\ to obtain (x = x 3 ) 



dr 



37 



Ax 



Araa 
9tt 



2-x {l-x)x 2(1-2;) 2 log(l-x) 



(2-xY 



{2-xf 



2 (1-x) log(l-a;) 



(16) 



which is the Ore-Powell result and is plotted in Fig. El 
Integrating Eq. (|16fl over x gives the total decay rate 



i^ — 



2 (tt 2 - 9) ma 6 
9^ 



(17) 



The low-energy photon spectrum results from the x — > 
limit of Eq. 



dr 



37 



d.r 



2ma b 
9tt 



I +<D{X2 < 



(18) 



which vanishes linearly with the energy of the radiated 
photon. 
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IV. EFFECTIVE THEORY COMPUTATION: 
MATCHING CONDITION 

The NREFT computation provides a systematic way 
of including bound state effects in the computation of the 
o-Ps decay amplitude. One computes the decay ampli- 
tude in QED. The NRQED Hamiltonian is constructed 
to reproduce the same amplitude, neglecting any binding 
effects. Once the NRQED Hamiltonian has been deter- 
mined, it can be used to compute the decay including 
binding corrections. 

The bound state dynamics is described by the Coulomb 
Hamiltonian for an e + e~ system in interaction with the 
quantized electromagnetic field: 

H = H o + Hmt 

TT - p2 a 
H — 

to r 

H int = -ii[*4, + * x ] B-ex E (19) 

with x = Xi — X2 the relative position of the pair of 
leptons, p! = p, p 2 = — p, and cr^, cr x the Pauli matrices 
acting on the electron and positron spinors (ji — e/2m). 

The interaction Hamiltonian H^t has been written in 
gauge invariant form in terms of the electric and magnetic 
fields. It can be obtained by a Goppert-Mayer transfor- 
mation of the usual Hamiltonian with a p — eA coupling 
(see e.g. Ref. [Tl|,[l2|)- The Goppert-Mayer transforma- 
tion removes the vector potential coupling in the covari- 
ant derivative, so that the electron and positron behave 
as neutral particles. The electromagnetic interactions are 
through multipole interactions with the electric and mag- 
netic fields. The electric and magnetic dipole interactions 
are shown in Eq. (|19l) . The higher multipoles arc higher 
order in v. Gauge invariance, and the computation using 
p • A interactions is discussed in Appendix. [5] 

The electric and magnetic fields E, B in H m t are evalu- 
ated in the dipole approximation, which is valid when the 
spatial extent of the positronium state, a = 2eto = 2/ma, 
is small with respect the wavelength of the radiated pho- 
ton (i.e. when k -C ma). The low energy photon is 
described by an ultrasoft field in NRQED, and the multi- 
pole ex pan sion follows from the NRQED velocity expan- 
sion |rlll3j. Our NRQED computation is valid provided 
the soft photon energy is small compared with ma. The 
NRQED result is valid even though the other two photons 
in the decay have energies of order to, because the anni- 
hilation into hard photons is a short distance coefficient 
in the NRQED Hamiltonian. It is possible to extend the 
NRQED computation to photon energies which are small 
compared to to but not ma by not using the multipole 
expansion. 

We shall use time-ordered (or "old-fashioned" ) pertur- 
bation theory (TOPT) in the effective theory calcula- 
tions, as it is more suitable for non-relativistic interac- 
tions. We recall that in TOPT vertices conserve only 
three-momenta and the virtual states are always on-shell. 
It is the violation of energy that characterizes the inter- 



mediate states rather than the off-shellness of the parti- 
cles, as in covariant perturbation theory. 

The Coulomb Hamiltonian Ho is the leading term in 
the velocity power counting. The kinetic energy and 
Coulomb potential are the same order in v. The en- 
ergies and wavefunctions of Hq are thus the Coulomb 
wavefunctions with reduced mass to/2. The electric and 
magnetic dipole interaction terms are treated as pertur- 
bations. The long distance part of the o-Ps decay ampli- 
tude involves these interactions. 

The p-Ps decay amplitude has no long-distance contri- 
bution, since both photons are hard, with energy to. The 
p-Ps annihilation amplitude is purely short-distance, and 
is given by computing the on-shell e + e~ annihilation am- 
plitude in QED as a power series in the lepton momentum 
p. For this paper, we need the first two terms in the v ex- 
pansion, which are given in Appendix. [S] Since the p-Ps 
decay amplitude is purely short distance, it can be writ- 
ten as a local operator in NRQED, with the hard photons 
treated as external sources from the point of view of the 
effective theory. The photons carry away the energy of 
order m produced in the annihilation, so the NRQED 
operator only depends on the scale ma. The first two 
terms in the v expansion computed in Appendix lAl give 
the annihilation operators in the Lagrangian (a sum on 
p is implicit) 

2 2 

(20) 

where ip p annihilates electrons with momentum p and 
x'-p annihilates positrons with momentum —p. We 
use the convention £0123 = +1 for the Levi-Civita ten- 
sor. The field strength tensors of the hard photons have 
been denoted by F to emphasize that they are external 
sources, and not dynamical fields in NRQED. The op- 
erators Eq. I|2(J|) allow one to compute the annihilation 
amplitude including all angular and polarization depen- 
dence, including that of the hard photons. One can in- 
stead integrate out the external sources F to get local 
four-fermion annihilation operators. These allow one to 
compute the angular and polarization dependence on the 
soft photon after averaging over all possible states of the 
hard photons. The four-fermion operators we will need 
are the S- and P-wave annihilation operators in Ref. pj, 

L = hr^l o^So) + /77(3 f o) 0( 3 P o ) + /77(3 f 2) 0( 3 P 2 ) 
m z m to 

(21) 

where the operators are [f| 

OCSo) = xlp^p |0) (0| 0t X _ p 

O( 3 P ) = ^ p (p-<x)0p|O)(O|0 P (p-<x)x-p 

0( 3 P 2 ) = X t - P (p ( V^)0 p |O)(O|0 P (p(V-'"))x-p 

(22) 



FIG. 3: Magnetic dipole graph for o-Ps annihilation. The 
solid square denotes the p-Ps annihilation vertex. 



FIG. 4: Electric dipole transition for o-Ps decay. The solid 
square denotes the p-Ps annihilation vertex. 



and the coefficients are given in Appendix A of Ref. Q 
Im/ r7 ( 1 5o) = 7ra 2 , 



Im/ 77 ( 3 P ) = 3vra 2 , 



Im/ 77 ( 3 P 2 ) 



4na 2 



(23) 



The annihilation Lagrangian can only be used to compute 
the square of the decay amplitude. 

The o-Ps decay amplitude has both short distance and 
long-distance contributions. The matching condition for 
the o-Ps decay amplitude in the effective theory is given 
by computing the difference between the QED o-Ps de- 
cay amplitude in the soft photon limit, Eq. JJJ, and the 
corresponding amplitude computed in the effective the- 
ory. The difference between the two computations gives 
the total short-distance o-Ps annihilation contribution in 
the effective Hamiltonian. We will compute the match- 
ing conditions using the annihilation amplitude Eq. <|2(J|) . 
and use it to compute the decay spectrum. In Sec. IVI Dl 
we show how the decay spectrum can be obtained using 
Eq. gl}. 



A. Magnetic Dipole 



The magnetic term in H lnt can induce a 1 3 Si —> 1 1 Sq 
transition between Coulomb e + e~ states. The effective 
theory graph for magnetic dipole emission is shown in 
Fig. and contributes to the long-distance o-Ps decay 
amplitude. The amplitude is: 



A.™ 



E 



i (0| \n) (n\ i fi [ar^ + a x ] ■ B |o-Ps) 

E Q — E n — E3 



. <0| A^ 1 1 1 fiT ) <1 1 J S r 1 1> [o-^ + <r x ] -B|o-Ps) 



E 



E p — E 3 



(24) 



where the only allowed intermediate state in the dipole 
approximation is the p-Ps ground state, with energy E p , 
so that E Q — E p = AEhfe. The energy of the emitted 
photon is E3 . Recoil effects in the energy propagator are 
suppressed by E 2 /m and are higher order in v. 

The magnetic dipole amplitude above is easily evalu- 
ated 



(l 1 S \i(J,[o-<f, + <T x ] -B|o-Ps) 



V2j 



-E 3 S 3 



(25) 

The quantity A^ 2 ^ is the e + e~ — > 27 amplitude with 
fermions at rest calculated in Appendix El see Eq. IA2I 
For e + e~ S'-wave states only the constant part of the 
A^ 2 ^ amplitude (Wo) contributes. When projected in 
the spin-singlet configuration (x^tp V%), it yields 



A (2 7 ) 



\/2i 



m 



[Si ■ e 2 + S 2 ■ ei] Vc(0). 



(26) 



To compute the o-Ps decay amplitude matching coef- 
ficient, one compares the effective theory result with the 
QED computation which neglects bound state effects. 
Neglecting the hyperfine mass difference in the energy 
denominator of Eq. (|25|l gives 



A, 



2m- 



tp o {0) X^3 ■ cr [81 ■ e 2 + S 2 ■ ei] ■ 



so we can define, in analogy with Eq. Q , 
V m = S 3 [Si ■ e 2 + S 2 ■ ei] . 
The limit of V m as P 3 — ► is 



(27) 
(28) 



V m = 



2 

+2 
+2 



ki ■ e 3 ) (e 2 ■ fc 3 ) - (ki ■ ^3) (^2 ■ e 3 ) 
( fe 3 ■ k^j (e 3 • e x ) - (k 3 • e x J (e 3 ■ ki 
ei ■ k 3 ) (e 2 ■ e 3 ) - (ei ■ e 3 ) (e 2 • k 3 



ei 



e 2 



ki. 

(29) 



B. Electric Dipole 

The electric dipole term — e x • E in H int can change 
orbital angular momentum by one unit, allowing for tran- 
sitions from ground state o-Ps to n 3 P 0)2 states (n ^ 1). 
We do not consider intermediate n 3 Pi states as they can- 
not decay to two photons ■ The effective theory graph 
for electric dipole emission is shown in Fig. ^ with am- 
plitude 



A, 



^ i (0| A^ \n) (n\ i ex • E |o-Ps) 



E — E n — E 3 



6 



and 

(n|iex-E|o-Ps) = e E 3 (n\ x • e 3 |o-Ps) . (31) 

For the matching computation, the energy separation 
between o-Ps and the intermediate states, E — E n ~ 
0(ma 2 ) is taken to be much smaller than the photon 
energy E3. In this case 



-4, 



n 

= -ie (o\ A (27) x • e 3 



(n\ x • e 3 |o-Ps) 



o-Ps 



(32) 



The two-photon amplitude A* 2,1 ' has terms which are a 
constant or linear in p. The constant term Wq yields the 
matrix element 

(O|x|o-Ps)=0 (33) 

since x?/> (x) vanishes at the origin. The linear terms 
Wi • p give 

(0|Wi-px-e 3 |o-Ps) 
which can be evaluated using 

-ftp ■ x + i (pV 



(34) 



p l x J 



(35) 



The o-Ps state is an S state, so to get a non-zero matrix 
element, the operator must have zero angular momen- 
tum. Only the first term contributes, and the p must act 
on x, otherwise the wavefunction vanishes at the origin. 
This gives 

(0|Wi-px-e 3 |o-Ps) = -i (0| Wi • e 3 |o-Ps) 

= - l ^(0) X t W 1 . e3 ^(36) 
Substituting the W term from Eq. (|A4|) we get 

o3 



A, 



so that 



ie 



ip o (0)x (e 3 -fci) (<r-fci) (ei-e 



+ (e 3 • 61) <t • e 2 + (e 3 ■ e 2 ) cr ■ ei 



(37) 



V e = 2 (e 3 • 61) e 2 + 2 (e 3 • e 2 ) e 1 + 2k 1 (e 3 ■ k^j (e 1 ■ e 2 ) . 



(38) 



The sum of the magnetic and electric dipole transi- 
tions in the effective theory, Eqs. I|29l) and l|38|) gives the 
full theory amplitude Eq. J7J). The matching condition, 
which is the difference of the two results, vanishes. Thus 
there is no additional three-photon annihilation term in 
the NRQED Hamiltonian, and the entire soft-photon o- 
Ps annihilation in the effective theory is due to electric 
and magnetic radiation followed by two-photon annihila- 
tion of p-Ps. This must happen, because the amplitude 
Eq. in the full theory is 0(1) as E 1 -> 0. Any local 
gauge invariant operator in the NRQED Hamiltonian de- 
pends on the soft-photon field-strength tensor, and gives 
an amplitude that is order 0{E~ I ). 



V. EFFECTIVE THEORY DECAY AMPLITUDE 
INCLUDING ENERGY DENOMINATORS 

The decay amplitude in the effective theory can now 
be computed without neglecting the energy difference be- 
tween o-Ps and the intermediate states. The matching 
calculation of the previous section shows that the entire 
contribution comes from the magnetic and electric graphs 
in Figs. EH 

The magnetic dipole amplitude including energy de- 
nominators is from Eq. I|24|) 



iA r , 



2m ; 



:ip o (0) 83 ■ cr [Si ■ e 2 + S 2 ■ ei] 



AE hfs - lu 
(39) 



while the electric dipole emission amplitude is 

( |A^)|n) (n|x- e3 |o-Ps) 
~ lAe = ^\ Eo-En-u, (40) 

where u> — E3 is the soft photon energy. The x • E op- 
erator connects to p-wave states, so only the Wi part of 
yj_(27) i s relevant, 



iA, 



E 



|Wi-p|rc) (nj x ■ e 3 |o-Ps) 

E„ — E n — LO 



(41) 



We can work out the matrix elements above with the 
help of the Wigner-Eckart theorem Define the re- 

duced matrix elements 

(n,l,0|z|o-Ps) = ((n,l|x|o-Ps)> 

(0|p z |n,l,0) = «0|p|n,l)), (42) 

where \nJtm) are the Ps states, and only 1=1 con- 
tributes. In terms of the reduced matrix elements, 
Eq. ijHJ) reads 



iA e 



E 



X f Wi • e 3 <j 



E — E n — lo 
x ((0\p\n, 1» ((ra,l|x|o-Ps)). 



(43) 



It is convenient to define 

to 



a m (uj) 
a e (ui) 

so that 
— iA m — 

-iA. = 



AEm s - u 



LO 



MO) 



E 



((O|p|n,l))((n,l|a;|o- Ps)) 

E n — E n — LO 



{44) 



2m 2 



V'o(0)a m x f S3 • cr [Si ■ e 2 + S 2 ■ ei] . 



e 3 • fci^ {cr ■ feij (ei • e 2 ) 



(e 3 ■ ei) cr ■ e 2 + (e 3 ■ e 2 ) cr ■ e 2 



0. (45) 
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From the equations above we already see that in the limit 
lu — > the effective theory amplitude A m +A e vanishes as 
lu, in agreement with Low's theorem. In the energy range 
Ai^hfs <C u> <C ma 2 the electric dipole amplitude still 
behaves as 0(lu), but the magnetic amplitude is 0(1). 
For the region above ma 2 both amplitudes are 0(1) and 
the Ore-Powell behavior is restored. 

Neglecting energy denominators, (i.e. taking the limit 
w 3> E n — E a ) the amplitudes a m and a e reduce to 



a m {uj — > oo) 
a e (uj — > oo) 



1 



^o(O) 



^((obKiH^i^io-Ps); 



(46) 



and we recover Eqs. (|?7jl and IpTTjl. 

To compute the spectrum from our effective theory 
calculation we first need to to sum over polarizations. A 
spin-1 state of o-Ps with polarization e is obtained by 
the replacement 



i 



e • er 



and the matrix element of cr is then 



(47) 



t). = -Ls-Trfa^] = V2e*. (48) 



The sum over polarizations is 

E 



131 



Sij k^kj, 



(49) 



for o-Ps and photons, respectively. Using these, it is 
straightforward to obtain the spin-averaged matrix ele- 
ment squared and summed over photon polarizations 



A, 



e e 

Jo 



l<Mo)l 2 {M 2 + kl 2 (2 + cos 2 0)} 



(50) 



(51) 



where cos 9 — k\ ■ k\. This agrees with Eq. II 1211 if one 
uses Eq. (gBJ. 

The phase space factor Eq. 1)15(1 can be written in terms 
of the variables cos# and x$ by using the relation (|13[) : 



X3 



l-X 3 



2 [l-^f (1-COS0)] 
^3 



■ dcos 9 dxq 



■ dcos 9 dx* 



(52) 



The above approximation is valid in the endpoint region 
where x% — > 0. The allowed region of integration for cos 9 
is — 1 < cos 9 < 1. 



As a etm are functions only of X3, the cos0 integration 
can be performed to give the differential rate (writing 
x 3 — > x): 



dT 

da; 



9tt 



l a m| 2 + - |a e | 2 



(53) 



The magnetic and electric amplitudes a m and a e are com- 
puted in the next section. 



VI. COMPUTING MATRIX ELEMENTS 

In this section, we compute the magnetic dipole am- 
plitude a m and the electric dipole amplitude a e . 



A. Magnetic Dipole 

For the magnetic term a m , the energy denominator can 
vanish, so we include the widths, 



LU 



where 



AEms — lo + ijo P 
7 



(54) 



A£hfs = —mo?, 

7o P = ^(r Q + r p ) w ima 5 . 

This gives the Breit-Wigner form 



(w - AE Ms f + 7o 2 p 
The magnetic amplitude behaves as 

f 1 to > AE hfs 



AE h 



u < AE 



(55) 



(56) 



(57) 



hfs- 



B. Electric Dipole 

We can rewrite the electric dipole amplitude, Eq. 



as 



^o(O) 



7?2 



E 



\n, 1, rn) (n, 1, to| 



o-Ps 



(58) 



where the term between brackets is the p-wave Coulomb 
Green's function G\ defined in AppendixlUland evaluated 
at energy —k 2 /m = E a — lu. In position space, using 

Eq. jnu, 



E 



(xjra, 1,to)(ti, Vrojy) 
+ fc 2 /m 



3(x-y) G 1 (x,y;k), 



(59) 
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gives the electric dipole amplitude 
Amu) 



a e {ui) 



MO) jo 



d«y*Gi(0,v;A)^ o (i/). (60) 



The integral in Eq. 1(60(1 can be evaluated numerically as 
a function of the soft photon energy u> using Eq. 1C14|) 
for i=l: 



Gi(0,y,k) 



mk 
"37 



e -% r(2- v) U(2- u,4,2ky), 



(61) 



where U is a confluent hypergeometric function, and 
v = am/ (2k). We note Eq. (|61|l for the p-wave Coulomb 
Green's function was used in Ref. [16(, with an additional 
factor of 6. Because of this discrepancy, we have given a 
derivation of Eq. (|HT|> in Appendix lO 
The electric amplitude behaves as 



a e (uj) 



1 



2cj 
mo- 



uj 3> ma 
uj <C ma 2 



(62) 



Note Added 



Recently, Voloshin 01 has shown that Eq. ((60(1 can be 
written in terms of a hypergeometric function, 



a e (uj) = 



(1-i/)(3 + 5i/) 
3(1 + i/) 2 

%v 2 {l-v) ( v-l 

+ ^7o \Tr, — \3 2-Fi 2 - z/, 1, 3 - i/; — — 

3(2 - i/)(l + j/) 13 \ z/ + l 



(63) 



C. Electric Dipole: Alternative Derivation 

The electric dipole amplitude a e can be obtained in a 
different representation by using the explicit form of the 
Coulomb wavefunctions. We follow the notation of Bethe 
and Salpeter and use the expressions found therein. 
Note that all quantities are dimensionless, with length, 
momentum and energy measured in units of l/(/xa), [ia 
and \xa 2 , respectively, where fi = m/2 is the reduced 
mass. 

For the discrete wavefunctions the dipole reduced ma- 
trix elements read 



;(n,l|z|o-Ps» = 



1 2 s n 7 (n- l) 2 ™" 5 
3 (n + l)2n+5 



1/2 



(64) 



IpIM» 



-#o(0) 



—(- 

4^3 \n 



5/2 r 



(n+l)(n-l) 



1/2 



(65) 



Define 



«0|p|n,l)) ((n,l|x|o-Ps)) _ 16 



n(n- I)™" 2 
(ri + 1)™+ 2 



(66) 

The continuum states Ryy(r) with energy W = k 2 /2 
are obtained from the discrete wavefunctions by the re- 
placement 



where 



n = T 



(67) 
(68) 



The exact relation reads 



Rni(r) -» (ik) 3/2 VT 



Rw(r), 



(69) 



with Z = 1 for e + e states. 

The completeness relation can be written as 



^ \n) (n\ + I dW\W) (W\ 



')■ 



(70) 



The sum over intermediate Coulomb states is given by 
the measure Eq. I|70|) . where the continuum states are 
normalized using the analytic continuation of the discrete 
states, Eqs. ((67169(1 . Note that the integral over n' in 
Eq. 1(70(1 has a non-trivial weight. 

The dipole reduced matrix elements for the continuum 
states read 

„ _ «0|pK,l))«n',lMo-Ps)) 



-#o(0) 

16 71/ r / i f\ f-\ 

= -TTT-, r^rr exp — ( 7T — 2 tan n)n\. 

3 (1 + n' 2 ) 2 ' 1 

(71) 

The p matrix element is given by analytic continuation 
of Eq. 1)65(1 . and the x matrix element by analytic con- 
tinuation of Eq. ((64(1 plus complex conjugation, since it 



involves 



rather than 



The electric dipole amplitude is (putting back the 
units): 



a e (u) 



E 



uj + E n 



dri 



o 1 - e- 2 ™' uj + E n > 



(72) 



with the energy differences from the ground state 



1 



E nn = -ma 1 r 



1 



1 



E n 'o = -ma 1 H — 



1 



(73) 
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Eq. iltjU|) using Eq. ijtjT|) for the Coulomb Green's function. 



D. Decay rate using local annihilation operators 

The decay spectrum Eq. I|53[l can also be computed 
using the annihilation Lagrangian Eq. i|21|) following the 
same steps as before. The magnetic dipole transition 
proceeds via S-wave annihilation, and the electric dipole 
transition via P-wave annihilation, and the diagrams 
have the form show in Fig. [SJ The decay spectrum is 



dr 

dx 



9tt 2 



Im/ 7T ( 1 5o) \a r 



1 



Im/ 77 ( 3 P ) \a e 



Im/ 77 ( 3 P 2 ) \a e 



(74) 



Using Eq. (|23J) reduces this to Eq. (£>3J). Of the 7 \a e \ 2 /3 
electric dipole contribution, the 3 Po contribution is |a e | 2 
and the 3 P2 contribution is 4 \ a e \ /3. 

In computing Eq. (|74|l . we have used the matrix ele- 
ment 



o-Ps 



xVVp |0) (0| 4>l^X-p 



o-Ps 



2 |^(o) | 2 44 

(75) 



0,01 



0. ' 







100 



E (Ryd) 



FIG. 6: Ratio of the o-Ps decay spectrum including bind- 
ing energy corrections to the Ore-Powell spectrum. The elec- 
tric dipole amplitude decouples for energies below a Rydberg, 
leaving only the magnetic contribution of 3/10. 



The magnetic contribution is given by the p-Ps resonance 
contribution. The photon decay spectrum in the p-Ps 
resonance region is shown in Fig. [7| At energies much 
smaller than the hyperfine splitting, the p-Ps state also 
decouples, and the decay rate vanishes as P 3 . 

The results are consistent with the Ore-Powell spec- 
trum and with Low's theorem. They include binding 
effects in a systematic expansion in powers of v. The 
modification of the decay spectrum also leads to a modi- 
fication of the total decay rate. Since the deviation of the 
spectrum is of order unity in a region of size the binding 
energy, and the decay spectrum is of order x, the total 
rate will be corrected at order a 4 by the effects we have 
computed. 



VII. RESULTS AND CONCLUSIONS 



The o-Ps decay spectrum is given by Eq. (|53[1 with 
\a m \ 2 given by Eq. lj5T)|) and a e given by Eq. ijrlUI) . The 
Ore-Powell spectrum is given by Eq. (|53[l with a m ^ e given 
by Eq. I|46|) . The ratio of the photon spectrum to the 
Ore-Powell value is shown in Fig. At energies large 
compared with the binding energy, a e and a m approach 
their values in Eq. 146(1 , and the spectrum approaches the 
Ore-Powell spectrum. The magnetic and electric dipole 
terms contribute in the ratio 3 : 7. Note that the ap- 
proach to the asymptotic value of unity is rather slow. 
At a photon energy of 100 Ryd, the ratio is 0.88. At en- 
ergies small compared with the binding energy, the elec- 
tric dipole transitions decouple, and one is left with the 
magnetic contribution of 3/10 of the Ore-Powell value. 
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E (eV) 

0.001 0.01 0.1 




0.001 
E (Ryd) 



FIG. 7: Ratio of the o-Ps decay spectrum including binding 
energy corrections to the Ore-Powell spectrum. The peak is 
the p-Ps resonance. The amplitude in this region is domi- 
nated by the magnetic dipole. For energies below the p-Ps 
resonance, the ratio to the Ore-Powell spectrum vanishes as 
E 2 . 



Pi 



lA/Wv k, 



P2 



FIG. 8: Two photon annihilation graph. 



APPENDIX A: PARAPOSITRONIUM DECAY 
AMPLITUDE 

The p-Ps-^ 27 decay amplitude is given by Fig. [8] 
We require the decay amplitude to first order in p. The 
incoming momenta are pi — (E, p), and p 2 = (E,—p) 

+ 0(p 2 ). The amplitude is 



with E = yj m 2 + p 2 



A (2l) 



-ie 2 v (p 2 ) L 



1 



h-h- 



(Al) 



We have not included any minus sign coming from Wick 
contractions for the initial e + e~ state. Including the 
crossed-graph, assuming the polarizations are transverse, 



and expanding in p gives 

A W = xHW + W 1 -p)<f> + O(p 2 ), 

where 



W a = — [Si -e 2 +S 2 - €1} 

2m 



w 1 = l 



(A2) 
(A3) 



a ■ feij (ei • e 2 ) fci 
(cr ■ e 2 ) ei + (a ■ e x ) e 2 



(A4) 



and we have used fci = — k 2 . 

APPENDIX B: GAUGE INVARIANCE 

One can compute in the effective theory using either 
x • E or p • A interactions. We have done the computa- 
tion using x • E. If one uses p ■ A interactions, then the 
gauge invariant momentum operator with our convention 
for the covariant derivative is p — eA. The p-Ps decay 
amplitude Eq. l|A2fl is then written in gauge invariant 
form as 

ffo + W r (p- eA) , (Bl) 
and the interaction Hamiltonian Eq. (|19|) has the term 

2e 



#int — 



-p A, 



(B2) 



instead of —ex • E. The electric dipole amplitude contri- 
bution from Eq. (|B2J| is 



2ei ^ <0| \n) i (n| p ■ e 3 |o-Ps) 



E 



Now 



E Q — E n — E3 



2i 



[H,x] = p, 

m 



(B3) 



(B4) 



so that the amplitude is 



oU (27) 



n > m x • e 3 



|o-Ps) 



E 3 {0\A (2 ^ \n) (n| x • e 3 |o-Ps) 



E n — E„ 



Eo 



(B5) 



The second term agrees with the dipole emission ampli- 
tude Eq. I|40|) using the x • E interaction. The first term 
gives the extra contribution 



ie(0 



A (2f) 



x ■ e 3 



o-Ps 



|Wi-e 3 |o-Ps).(B6) 



In addition, one has the amplitude from the — eA term 
in Eq. 

- e (0| Wi • e 3 |o-Ps) (B7) 

which exactly cancels Eq. (|B6ll . so the net contribution 
is gauge invariant. 
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APPENDIX C: COULOMB GREEN'S FUNCTION to get 



The Coulomb Green's function satisfies the differential 
equation 



H„ + — )G(x,y,fe)=<5(x-y), (CI) 
m 



with the Coulomb Hamiltonian for positronium 



m r 

In terms of a complete set of states in position space 

i) n (x) V* (y) 



(C2) 



G(x, y ,fc) = Y, 



E n + k 2 /m 



The Green's function admits a partial wave decomposi- 
tion: 



G(x,y,fe) = ^(2£+l)(^) £ P,(x-y/xy)Ga^2;,fc) 

(C4) 



1=0 



where Pg(x) are the Legendre pol yno mials. The partial 
waves of the Green's function are 



Gi (x,y, k) 

= mfc (2k) u V Lf+\2kx)Lf+\2ky)r\ 

2tt J ^ (r + £+l-j/)(r + 2£+l)! 

r— 

(C5) 

with the parameter 

am 



and the associated Laguerre polynomials defined by 



(C6) 



L k r (x) = 



° x -~ k ' d 



-\ V Ax 



e - x x r+k . 



(C7) 



The Green's function partial waves have a more com- 
pact expression when one of the arguments is zero |lfi| . 
Set y = and use 



(r + k)\ 
r\k\ 



(C8) 



G e {x, 0, k) 

j OO 

— (2k? 1 e~ kx V 

2tt K ' ^ (r + l+l-v)(2l+l)V 

r— 



Lf +1 {2kx) 



(C9) 



The generating function for the associated Legendre 
Polynomials [2(| is 



(C3) Also 



OO 

r=0 



■U5/(l — S i 



1 - s) a+1 



1 



ds 



r+b-1 



r + b jo 

The identities (|C10|) and I|C11|I allow us to write 



(CIO) 



(Cll) 



^ r + b 

r— 



-us/ (1 — s) 



Changing variables to £ = s/(l — s) gives 



E 

r=0 



Lt(u) 



a—b^—ut 



dtt b -\l + t) a - b e 



T(b)U(b,a+l,u), 



(C13) 



where we have adopted the usual definition of the conflu- 
ent hypergeometric function U(a,b,z), see e.g. Ref. [21j . 
Applied to the sum in Eq. EH), the result (|CT3|) yields 

G t (x,0,k) 

mk 2l _ kx T(£ + 1 - v)XJ{l + 1 - ^, 21 + 2, 2kx) 
~ 2n {2k) 6 (2* + 1)1 

(C14) 

which differs from Ref. 

[3 by the (2^+1)! factor in the 

denominator. 
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